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T REFACE 


The investigation of frontogenesis in a three- 
dimensional linear velocity field was commenced in the 
fall of 1947 at the U. S. Haval Postgraduate School in 
partial fulfillment of the requirements for the degree 
of faster of Science* 

The two-dimensional development of frontogenesis 
by Petterssen (5*6) aroused the desire to visualize and 
to attempt to portray what would occur in an extension 
of his definition of frontogeneBis into three dimensions. 
This paper covers the main essentials in the development. 
Time limitations have prevented a more exhaustive 
investigation. 

The author wishes to express his appreciation for 
general guidance to Professor W. D. Duthie, and for 
assistance by the staffs of the departments of Aerological 
Engineering and Mathematics of the U. S. Haval Postgraduate 
School. 



TABLE OF CONTESTS 


Page 


CERTIFICATE OF APPROVAL i 
PREFACE ii 
LIST OF ILLUSTRATIONS iv 
TABLE OF SYMBOLS AND ABBREVIATIONS V 
CHAPTER 

' I. INTRODUCTION 1 


II. LINEAR FIELDS OF MOTION 

1• Linear Vector Function of Position 
2. The Equations of Linear Fields of 
Motion 

3# Streamlines 
4. Translation Field 
5* Divergence Field 

6. Deformation Field 

7. Rotation Field 

8* The General Field 



III. FRONTOGEHESIS II GENERAL 15 

9. Definition of Frontogenesis 15 

10* Frontogenesis in Linear Fields of 

Motion 15 

11. Frontogenetical and Frontolytical 

Sectors 16 


IV. FRONTOGENESIS IN ATMOSPHERIC WIND SYSTEMS 21 


12. Preliminary Discussion 21 

13. Frontogenesis in the Vicinity of 

Cols 22 

14. Frontogenesis in the Vicinity of 

Cyclones and Anticyclones 25 

15* Frontogenesis in the Vicinity of 

Troughs and Wedges 28 


V. CONCLUSIONS 30 

BIBLIOGRAPHY 32 

APPENDICES 

I. 1. Linear Vector Function of Position 33 

2. The Equations of Linear Fields of 

Motion 33 

II. Solutions of Linear Differential 

Equations with Constant Coefficients 37 

III, Frontogenesis in Linear Fields of 

Motion ai 


(111) 



LIST 03? ILLUSTRATIONS 


Page 

1* Deformation field with Convection. 8 

2, Deformation field with Subsidence, 9 

3, Spheroids of Deformation. 10 

4, General field for the Heal Case* 12 

5* General field for the Complex Case. 14 

6. Frontogenetical and Frontolytical Sectors 

for C Positive. 17 

7* frontogenetical and Frontolytical Sectors 

for C* Negative. 19 


(It) 



TABLE OF SYMBOLS AND ABBREVIATIONS 


V 

i,j»k 

r 

x,y,z 

u»v»w 

u 0 ,v o ,w* 

u x ,(u 3 ), etc. 

§>(§c) 

®tC. 

K lf (Kj, etc. 

P, * Pz » ft 

A, t-0-2*-^-3 

a+bi, etc. 

In a 

l»m,n 

f 


vector velocity 

unit vectors along the x,y»z axes, 
respectively, in a right handed coordinate 
system 

position vector 

rectangular coordinates 

scalar components of vector velocity, V, 
along the x,y»z axes, respectively 

constant velocity components at the point 

( 0 , 0 , 0 ) 

partial differentiation of u with respect 
to x,(y), etc. 

tensor, (conjugate of the tensor) 

angular velocity of rotation about the 
x, (y), etc. axis 

arbitrary constant arising from the 
process of integration 

arbitrary constant multipliers 

roots (real and/or complex) of a cubic 
equation 

linear functions of x,y,z 

complex number of real part, a, and 
imaginary part, b 

natural logarithm of a (to the base e) 
direction cosines 
frontogenetical function 


S 


any conservative property that may be 
observed in the free atmosphere 



absolute magnitude of the gradient of 
property* 




frontogenesis. 

frontolysis. 


Fg 

FI 



I, IBTROBUCTIOH 


The first available source of information on 
frontogenesis is by Bergeron (1). He gives a limited 
treatment of the deformation field only, in three 
dimensions. Bjerknes (2) later treats the case of the 
deformation field in two dimensions. The first full 
treatment of frontogenesis is given by Petterssen (5>,6). 

His treatment is limited to two dimensions, however. 

This paper is a generalization of the two- 
dimensional theory of frontogenesis in a linear velocity 
field as developed by Petterssen (5). 

The most general homogeneous linear velocity field 

in three dimensions contains nine constants to describe 

the motions of translation, divergence, deformation, and 

rotation. The two-dimensional specialization of the 

general case contains five constants to describe the same 

motions. ; 

* 

An investigation of the possible frontogenftical and 
frontolytical sectors determined by the surface for which 
the frontogenetical function, f, is equal to zero yields 
the following! 

(1) The surface for which f is equal to zero, if it 
exists, is an elliptic cone with vertex at the origin of 
the coordinate system. 

(2) The wind field, as described by the nine constants 
mentioned above, determines whether or not the surface for 



which f equals zero will exist* If it does exist, certain 
of these constants determine just what the shape, size, 
and orientation of the conical sectors will be. 

(3) If the direction of tne gradient of property 
being considered lies (a) within the frontogenetical 
sector, (b) on the surface for which f equals zero, or 
(c) outside the frontogenetical sector, there will take 
place (a) frontogenesis, (b) no change, or (c) frontolysis, 
respectively. 

The conservative property ttsed in this paper is potential 
temperature whose equipotential surfaces are nearly horizontal 
and whose gradient is nearly vertical. Instability can 
cause the gradient of property to depart from the vertical. 

An application of the theory of frontogenesis is made 
in each of the following types of wind systems observed in 
the free atmospheres cols, cyclones, anticyclones, troughs, 
and wedges. In each of these cases only the normal situation 
is treated in this paper. It is found that in the case of 
a col with convection at its center frontolysis will occur 
at the center, while in a col with subsidence at its center 
frontogenesis will occur at the center. In the case of a 
cyclone with convection at its center and attendant inflow 
in the horizontal plane frontolysis will take place at the 
center. In an anticyclone with subsidence at its center 
and attendant outflow in the horizontal plane frontogenesis 
will take place at the center. Troughs and wedges in this 
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paper are considered as special cases of cyclones and 
anticyclones, respectively* JTontogenesis or frontolysis 
occurs as is indicated for the cyclones or anticyclones, 
if the direction of the gradient of property is assumed 
nearly vertical. 

Variations in the frontogenetical effect can occur 
in the above svstems In the area around the centers of 
the systems where marked convection or subsidence can 
cause the direction of the gradient of property to 
depart from the vertical. 


(3) 



II. LISE JUS FIELDS OF HOTIOS. 

1* Linear Vector Function of Position. 

For the velocity# V, the linear vector funotion of 
position can he written as 

( 1 - 1 } V = \4 + /t • v V. 

In Cartesian coordinates this becomes 

us U. + X 0 * + y U 3 + z o 4 

/ (1-2)) crscr® + x ir» + yu^+Zu** 
or* ur. + X ur* + tj ur 3 +2 ur fe . 

This development follows that of Prandtl and Tietjene 
(7) and is given in section 1 of ippendix I. 

The equations above define the velocity field in the 
vicinity of a point where the second order and higher 
terms of the Teylor series are negligible as compared to 
the linear terms. 

2* The Equations of Linear Fields of Motion. 

The equations (1-2) are transformed in section 2 of 
Appendix I so that the equations of linear fields of motion 
are 

0 s Oo+Ax + Dy + Eg 
(2-1) U , “U’o+Bj-DX + / r 2 

ur * ur 0 + C 2 -£x-/ : ’'jj 

or, rewriting, 

Oo + (A+ B+c) x -(B + Ox + Dy +£* 

(2-2) cr* (Jo ■+ (A + B+- c) j “(A+C)«j - Dx +A 2 

ur= of. +-CA + B-*-C)2“(A^ B)Z - £ X - Acj , 



/ 
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where; u 0 ,v 0 ,w B are the constants of translation; (A+B+C) 
is the constant of divergence; (B+C), (A+C)» (A+B) are the 
constants of deformation; D,E»F are the constants of 
rotation; end A,B,C are the constants of extension along 
the x,y,z axes* respectively* 

3* Streamlines* 

A streamline is a line for which the tangent at any 
point gives the direction of the fluid motion at that 
point* 

The equations of section 2 will be discussed in terras 
of the streamlines for convenience. The differential 
equations of the streamlines are 


(3-1) = . 

O w or 


These equations will be used to study, first, the separate 
components of translation, divergence, deformation, and 
rotation, end then the general case in which all of the 
components are included* 

4* Translation Field. 

Considering only the component of translation of 
equations (2-2), solutions of the differential equations 
give 


(4-1) 


x- k. , 3- *z. _ a -*3 

lT 0 uJo 


where u 0 *v. »w 0 ere the direction numbers of a family of 
straight, parallel streamlines through the points 
determined by the arbitrary constants K, »X 3 * This is 
a rectilinear field of motion. 

(5) 



5* Divergence Field, 

Considering the component of divergence only* 
solutions of the differential equations give 
(5-1) K= Kj-y = M . 

The streamlines are straight, radial lines emanating 
from the origin of the coordinate system* Mathematically, 
the origin is a source of fluid if the divergence is 
positive, and it is a sink for the fluid if the divergence 
is negative (convergence)* These cases are (j+B+C) 
greater than* or less than, zero, respectively* 

In this paper, divergence refers to mass divergence 
which can express itself only as a change in the fluid 
density. In the two-dimensional case as presented hy 
Stewart (8), the divergence is a "velocity" divergence 
where the third dimension acts as the source or sink of 
fluid* and where the mass divergence is equal to zero. 
Such a velocity divergence is represented as the three- 
dimensional deformation field* 

6* Deformation Field. 

If only the component of deformation is considered, 
solutions of the differential equations are 

X«j2.= K n 

(6-1) 



where -(e+e) = A 

- (A + c) = B 

- (A + 8) = c 


l 
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since A+B+C equals z&io. Choose axes so that A is always 
positive, B is always negative, and C is either positive 
or negative* Let this sign convention he denoted as 

As +A ‘ 

(6-3) B = - 0' 

c — ± c' . 

This convention causes the x axis to he an axis of 
dilatation, the y axis to he an axis of contraction, and 
the 2 axis to he an axis of dilatation if 0* is positive 
and sn axis of contraction if G* is negative* When A* 
and C # are positive with B* negative, this causes the xz 
plane to he a plane of outflow and the y axis an axis of 
inflow* When is positive with B* and G* negative, 
this causes the yz plane to he a plane of inflow and the 
x axis an axis of outflow* In applying the above solutions, 
or the general solutions that follow, to actual examples, 
it is necessary to follow this sign convention established 
in order that succeeding developments will apply. 

From the above solutions it is seen that two stream¬ 
line patterns can be obtained depending on the sign of C, 

If C f is positive (negative) the streamline pattern will 
be similar to that shown in figure 1 (2)* Positive 
(negative) C* represents outflow (inflow) along the z 
axis, with outflow along the x axis and inflow along the 
y axis* For C # positive (negative), this i© the 
deformation a sphere would undergo in ohanging into an 
oblate (prolate) spheroid. Refer to figure 3* 
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FtGURE 1. 



Figure 2 . 
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Fiaoftc 3. 
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7. Rotation Field. 

Considering only the component of rotation, solutions 
of the differential equations are 

(7-1) X 2 + ^ +2*=K to 2 + 

uhere K.ois greater than K„ for streamlines to exist. A 
simultaneous solution of these equations gives a family 
of circular streamlines all with a common axis of rotation 
defined by the direction numbers F ,-£,D. 

8* The General Field. 

In order to solve the differential equations when 
all the components, except translation, are considered it 
is necessary to make the equations exact* There exist 
multipliers such that the numerators are constant multiples 
of the differentials of their respective denominators. 

These multiples are designated p f , p^ » p a « They are the 
roots of a cubic equation, and as Buch there may he all 
real roots or one real root and two conjugate complex 
roots. This derivation is due to Ford (3) and will he 
found in Appendix II* 

For the case of all real roots the solution is 

(8-D 

where xi,,xi z ,n 5 are linear functions of x,y»z with real 
terms only. A simultaneous solution of these equations 
gives skew space curves for the streamlines similar to 
those shown in figure 4. The sign convention of section 
6 applies here to the oblique planes and axes. For A* 
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and C* positive with B* negative, the xz plane is a plane 
of dilatation; the x*z* plane is a plane of outflow; the 
y axis is an axis of contraction; and the y* axis is an 
axis of inflow. For A* positive with B* and C* negative, 
the yz plane is a plane of contraction; the y^* plane is 
a plane of inflow; the x axis is an axis of dilatation; 
and the x* axis is an axis of outflow. 

For the case of imaginary roots, it is shown in the 
Appendix that the solution is of the form 


( 8 - 2 ) 


% inS 1,= in KU+ '/z K [(cx+e^+gz)* + (d x + -fv hz)] 

fl 


f, cx+ey + <3* 


where il,, is a linear function of x,y,z with real terms 
only, A simultaneous solution of these equations gives 
skew space curves of a general helical, or spiral, shape* 
Refer to figure !?* 
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III. FRONT0GEN3SIS IN GENERAL. 

9. Definition of Frontogenesis* 

The definition of frontogenesis as given by 
Petterssen (5) will be used. Essentially, it is the 
tise rate of change of the gradient of a conservative 
property which can be written 

f s b s| . 

Petterssen has developed this definition further, 
so that it becomes 

(9-2) f =-vs-(vs- 7)V / Y 7s| 

10. Frontogenesis in linear Fields of Motion. 

When equation {9-2} is applied to the equations 
for linear fields of motion (2-1), the following 
equations result: 

(10-1) -f = - \*S l[A ft + 6 v^-t-Cn 2 ] 

(10-2) f = |vS|[(6 + C)F+(A4c)>n Z + (A + B) n"- (A + B + c)] 

where l,m,n are the direction cosines of the gradient 
of property. The intermediate steps in the derivation 
are given in Appendix III. These equations show that 
the frontogenetical effect is independent of both 
translation and rotation. Divergence deducts from the 
frontogenetical effect} that is, it produces frontolysis. 
Convergence adds to the frontogenetical effect} produces 
frontogenesis. The bearing deformation has on the 
frontogenetical effect depends on the values of the 
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deformation constants. It is obviouB that if there 
were no velocity of extension occurring in the atmospheric 
wind field* there would be no frontogenetical effect. 

The Burface in which the absolute magnitude of the 
gradient of property is a maximum coincides with the 
surface of frontogenesis. 

11. Frontogenetical and Frontolytical Sectors. 

The equations (10-1) and (10-2) show that the sign 
of f is determined by the quantity within the brackets. 

An investigation of the surface for which the frontogenetical 
function is equal to zero yields a description of the 
frontogenetical (f greater than zero) and frontolytical 
(f less than zero) sectors with relation to the gradient 
of property. 

For f equal to zero* 

( 1 1 — 1 } A £%• 6 >n* +• C tx L =• O ^ 

or • 

(11-2) (B + C) £*+ (A+C)-m+(A+ 8 )^-(a+B + C) =o 

These equations represent an elliptic cone with vertex 
at the origin of the coordinate system. Prom the sign 
convention established in equation (6-3)» two orientations 
of the cone are possible. If C* is positive, the y axis 
(axis of contraction) will be the axis of the cone. The 
frontogenetical sector will be contained within the cone* 
end the frontolytical sector will be all the space outside 
the cone* as shown in figure 6. If G* is negative* the 
x axis (axis of dilatation) will be the axis of the cone. 

The frontolytical sector will be contained within the cone, 
and the frontogenetical sector will be all the space 




outside the cone# as shown in figure 7* 

At this point it is to he emphasized that the cones 
are symmetric with respect to the axes of dilatation or 
contraction and not with respect to the axes of outflow 
or inflow. This holds from the previous developments 
where it was shown that rotation separates the axes of 
inflow and outflow from the axes of contraction and 
dilatation* respectively# while rotation has no effect 
on the orientation of the frontogenetical and frontolytical 
sectors. 

If divergence predominates over the deformation# 
there will he no sectors# and no frontogenesis (all 
frontolysis) will occur. However* if convergence 
(negative divergence) predominates over the deformation# 
again# there will he no sectors# hut frontogenesis will 
oocur everywhere within the region being considered if 
the region is linear and if the space rate of change of 
the orientation of the gradient of property is negligible. 

Por each of the chosen values of A#B#C there are 
eight possible orientations of the gradient of property, 
since l,m,n are each squared terms in the frontogenetical 
function. 

If deformation predominates over divergence# or 
convergence, the orientation of the gradient of property 
determines whether or not frontogenesis will occur. If 
the direction of the gradient of property lies (a) within 
the frontogenetical sector# (b) on the surface for which 


(18) 




Figure 7 . 
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f is equal to zero, or (c) outside the frontogenetical 

(within the frontolytical) sector, there will take place 

(a) frontogenesis, (b) no change in the frontogenetical 

effect, or (c) frontolysis, respectively. These 

considerations may be applied in any of the eight octants 

of the coordinate system, since a portion of the sectors 

is present in each octant, and the direction of the 

gradient of property may appear in any octant* It is to 

be noted that these sectors, as defined, are interchanged 

with those as given by Petterssen (5*6)* His sectors 

are related to the isolines of the conservative property, 

while the sectors herein are related to the gradient of 

♦ 

the property. 

If potential temperature is used as the conservative 
property in the frontogenetical function, the direction 
of its gradient is nearly vertical. However, under 
conditions of marked subsidence or convection, the 
direction of the gradient would depart from the vertical. 
This point is brought out in the chapter on atmospheric 
wind systems. 
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IV. PROHTOGSEBSIS IS ATMOSPHERIC WIED SYSTEMS. 

12. Preliminary Discussion. 

The dominant types of atmospheric wind systems 
observed in the free atmosphere are cols* cyclones, 
anticyclones, troughs, and wedges. In order that the 
mathematical developments may apply to these systems, 
one of the coordinate planes will be assumed nearly 
horizontal, and then the coordinate axis normal to this 
plane will be nearly vertical. 

The original derivation of the equations of linear 
fields of motion in section 1 required that the discussion 
be limited to the immediate vicinity of a point so that 
the second order and higher terms of the Taylor expansion 
would be negligible in comparison with the linear terms. 

In the discussion to follow it must be assumed that the 
entire region under consideration is a linear velocity 
field wherein the velocity components are linearly 
dependent on the distance from the point in question. 

This point is usually the center of the atmospheric system 
under consideration. 

In section the possible interpretations of 
divergence were discussed mathematically. There are no 
pressure fluctuations observed in the atmosphere that 
represent pure mass divergence or convergence in stable 
systems of the size being considered. The discussion 
to follow is limited to stable systems and therefore 
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will assume that divergence is negligible in comparison 
to deformation and rotation in the shaping of the wind 
systems observed, 3y stable systems are meant those 
systems that retain their identifying features for a 
reasonable length of time. For these assumptions* the 
frontogenetical function can be simplified to 
(12-1) f- \*S\(-At 2 -d^-cW-) , 

since A+B*C equals zero. Then* in atmospheric systems, 
choose axes so that A is positive, 3 is negative, and 
so that C is positive or negative (equation ( 6 - 3 )), 

The conservative property to be used in this paper 
is potential temperature. Since the equipotential 
surfaces are nearly horizontal, the gradient of potential 
temperature will be nearly vertical. As one approaches 
the center of an atmospheric system with marked convection 
or subsidence, the gradient will depart from the vertical 
by an amount depending primarily on the constant of 
extension in the vertical. At the very center, however, 
the direction of the gradient must return to the nearly 
vertical direction, 

13, Frontogenesis in the Vicinity of Cols. 

A col is a region of relative calm surrounded by 
four distinct systems{ two diametrically opposite cyclonic, 
low pressure, systems and two diametrically opposite 
anticyclonic, high pressure, systems, Petterssen ( 6 ) 
distinguishes three types of oolsi (a) without vorticity, 
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in which the relative sizes of the high and low systems 
with respect to the col are the seme* (b) with anti- 
cyclonic vortioity in which the relative sizes of the low 
By8terns predominate over the relative sizes of the high 
systems, and (c) with cyclonic vorticity in which the 
relative sizes of the high systems predominate over the 
relative sizes of the low systems* 

In a col the constants of extension predominate over 
the constants of rotation* The three-dimensional case 
superimposes either ascending or descending motion on 
the accepted two-dimensional case as given by Petterssen 
(5*6) and Stewart (8}« 

For the xy plane to be nearly horizontal 
(13-1) |A + 0l»|£| and lA-*-6|»|fl 

or | C | » | £ | and I c'l » I ? I . 

And for the rotation to be a minimum about the z axis 
(13-2) |D|«|8 + C| and (Dl«lA+c| 

or |D|«|A'| and 10|«I &'I . 

When these conditions are met* then two types of cols 
are possible. If C* is positive there is ascending motion, 
or convection* at the center, while if C* is negative 
there is descending motion* or subsidence, at the center. 
The above restrictions cause the equations for linear 

fields of motion (2-2) to take the form 

O - Go + a' * 

(13“3) If =T Ui - 5 / *j 

ur= of, ± C'z , 

( 23 ). 



if tJie smaller magnitude terms are considered negligible. 
This is a special case of the general field. Essentially, 
these are the equations used to obtain the solutions of 
section 6 for a pure deformation field. Slight amounts 
of rotation vould cause the coordinate system to become 
oblique, as shown in figure 4. 

It must be recalled from section 10 that the 
frontogenetical sectors are determined only from the 
constants of extension. Since, in the case of a col, 
the constants of extension are large compared to those of 
rotation* figure 1 superimposed on figure 6 represents 
frontogenesis in a col with convection* end figure 2 
superimposed on figure 7 represents frontogenesis in a 
col with subsidence* 

Figure 1 on 6 shows that the frontogenetical cone 
is symmetrical with respect to the axis of contraction 
and that the frontogenetical sector is contained within 
the cone. It is to be recalled from section 11 that the 
cone is not affected by rotation even though the axis of 
inflow does move away from the axis of contraction* and 
even though the plane of outflow does move away from the 
plane of dilatation* when rotation is present. Since the 
direction of the gradient of property is nearly vertical 
at the center* frontogenesis would not take place in this 
case with convection. From the preliminary discussion 
of the preceding section, in the area around the center 
of the col the direction of the gradient of property may 
depart sufficiently from the vertical for it to lie within 
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the frontogenetical cone in which case frontogenesis 
could occur in the direction of the axis of contraction* 

Figure 2 on 7 shows that the frontogenetical cone 

is symmetrical with respect to the axis of dilatation 

and that the frontogenetical sector is all the space 

outside the cone. Again# since the direction of the 

gradient of property is nearly vertical, frontogenesis 

could take place quite readily in this case* However, 

around the center of the col the direction of the gradient 

of property may depart sufficiently from the vertical so 

that it would lie within the frontolytical sector* 

14. Frontogenesis in the Vicinity of Cyclones and 
Anticyclones. 

A cyclone, as used herein, refers to a northern 
hemisphere extretropical cyclone shout which the motion 
is circular, or nearly circular, in a counterclockwise 
direction about a center of low pressure. An anticyclone 
refers to a northern hemisphere region of relatively high 
pressure about which the motion is circular# or nearly 
circular# in a clockwise direction* 

A cyclone requires that the constants of extension 
predominate over the constants of rotation in two of the 
coordinate directions while the reverse must be true in 
the third direction* The normal cyclone has inflow in 
the horizontal plane with attendant convection (or outflow) 
in the vertical direction* The sign convention established 
in section 6 requires that the yz plane be nearly 
horizontal and that the x axis be nearly vertical. In the 


{ 25 ) 



notation being used, these conditions are 

(14-1) I A' | » ID | l | A / 1 >•> 1 £. | 

which causes the yz plane to be nearly horizontal, and 

(14-2) |f r |»| d'l $ \f\ » l C*\ 

which causes the rotation about the x axis to be a 
maximum. Further, for cyclonic (counterclockwise) rotation, 
F must be negative. 

The above restrictions cause the equations for linear 
fields of motion (2-2) to take the form 

O = Oo + fit x 

(14-3) LS=: CTo - ^2 

US - (jSo -h 

if the smaller magnitude terms are considered negligible. 
This is a special case of the general field. 

"Figure 5 superimposed on figure 7 shows that the 
frontogenetical cone is symmetrical with respect to the 
axis of dilatation (also axis of outflow in this case) 
and that the frontogenetical sector is all the space 
outside the cone. Since the direction of the gradient 
of property is nearly vertical, frontogenesis would not 
occur in this case. It is observed in the free atmosphere 
that as one approaches the center of a cyclone, frontal 
discontinuities are difficult to locate. Again, from 
the preliminary discussion of section 12, in the area 
around the center of a cyclone the direction of the 
gradient of property may depart sufficiently from the 
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vertical for it to lie within the frontogenetical sector, 
in which caBe front ogenesis could occur. 

An anticyclone requires that the constants of 
extension and rotation have the same relative magnitudes 
as they do for a cyclone. However, for an anticyclone 

V 

the direction of the rotation is reversed* A normal 
anticyclone has outflow in the horizontal plane with 
attendant inflow (or subsidence) in the vertical direction. 
The sign convention of section 6 requires that the xz 
plane be nearly horizontal and that the y axis be nearly 
vertical. These conditions are 
(14-4) |B'|»|D| & lB'J»lF r l 

which causes the xz plane to be nearly horizontal, and 

(14-5) UI»|A'I & |£|»|c'| 

which causes the rotation about the y axis to be a 
maximum. Further, for anticyclonie rotation, E must be 
negative. 

The above restrictions cause the equations for linear 
fields of motion (2-2) to take the form 

O « Oo - 

(14-6) U'=cr.-6 / y 

Uf = US 0 + £ % 

if the smaller magnitude terms are considered negligible* 
This, too, is a special case of the general field. 

The frontogenetical cone is symmetrical with respect 
to the axis of contraction (also axis of inflow in this 
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case), and the frontogenetical sector is contained within 
the cone. Again, since the direction of the gradient 
of property is nearly vertical, frontogenesis could occur 
c^uite readily in this case. From previous discussions 
it can be seen that in the area around the center of 
an anticyclone frontolysis could take place. 

15* Frontogenesis in the Vicinity of Troughs and WedgeB. 

A trough is an elongated area of relatively low 
pressure extending from the center of a cyclone, while 
a wedge is an elongated area of relatively high pressure 
extending from the center of an anticyclone. 

The essential difference between troughs and wedges 
and cyclones and anticyclones is that the region under 
consideration is no longer the center of the system but 
is an area outside the influence of the center itself. 

The frontogenetical cone has the same orientation 
it had in the case of the cyclone, since translation 
has no effect on the frontogenetical function. As in 
the cyclonic case, the direction of the gradient of 
property with respect to the frontogenetical sector in 
a trough would cause frontolysis to occur, assuming the 
direction to be nearly vertical. 

A wedge is analogous to an anticyclone and hence 
is conducive to frontogenesis. 

Until further investigation can be made into the 
actual magnitudes of the direction cosines of the gradient 
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of property with respect to the linear velocity field 
constants in troughs and wedges, the following must he 
inferred from the discussion of cyclones and anticyclones. 
Since troughs and wedges are being considered as extended 
areas outside the vicinity of the centers of cyclones 
and anticyclones, respectively, it is quite possible 
that the direction of the gradient of property is no 
longer vertical so that frontogenesis can occur in a 
trough and frontolysis can occur in a wedge* For example, 
figure 6 shows that frontogenesis will occur more readily 
in the direction of the y axis than the z axis due to 
the elliptical shape of the conical sector* 

Fetterseen (6), pg. 258, gives an informative 
description of the relationship between the movements 
of systems and the movements of fronts. In circular 
or hyperbolic systems it is necessary for the system 
to move with nearly the speed of the front for the front 
to maintain itself. Should the system be stationary, the 
front would be subject to first frontogenesis, then 
frontolysis, and then frontogenesis again, etc. 
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V. COUGLUSIOirS 


Prom the foregoing developments it can "be seen that 
this paper has no more than probed the surface of the 
topic of Frontogenesis in a Three-Dimensional Linear 
Velocity Field. In the time available it has not been / 
possible even to parallel all of the topics presented 
by Petterssen (5*6) for the two-dimensional case. It 
is the intention of the author to carry forward the 
investigation commenced herein* not only from a purely 
theoretical viewpoint* but also along practical lines 
to see whether or not another useful forecasting tool 
can be devised# 

At its present stage of development this paper should 
give the reader some idea of the processes that can occur* 
or are likely to occur, in the formation* persistence* or 
dissipation of fronts in the free atmosphere. Petterssen 
(6) gives a complete description of the general fronto- 
genetical regions of the world. 

Further investigations should include: (a) research 
toward finding a more conservative property than potential 
temperature, (b) the percentage of the frontogenetical 
sector as related to the total solid angle about the 
origin of the coordinate system for different values of 
the constants of extension, (c) the large changes in the 
direction of the gradient of property in the vicinity 
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of unstable convective cells and the effect on frontogenesis, 
(d) the effect of variations in the relative magnitudes 
of A,B»C, (e) applications to actual synoptic situations. 
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APHBHDIX I. 


1* Linear Vector Function of Position* 

Let a fixed point in space have the coordinates 
(0,0,0) in a rectangular coordinate system (x,y,z)# 

And let the velocity at the fixed point he 
(A 1-1 ) Vo a t Uo + j CTo + KuTa 

where i»j»k are the unit vectors in the x,y,z directions 
respectively# Then the vector velocity, V* in the vicinity 
of the fixed point can he represented as a vector function 
of position hy an expansion of V in a Taylor series: 

(A 1-2) V = V« + A- * v V + ~r * v 2 V h- 

where v (del) has the standard mathematical definition 
and where the position vector, r, is 
(A 1-3) A- = tX + j <j + K2 . 

In Cartesian coordinates, neglecting second order and 
higher terms, V can he written as 

U * Uo + X + <j U^J +• £ o* 

(A 1—4-) W = cTo + X Wk + 'jCT^ + H.Ua 

(jlT= u/o + X ur* + ^ u£, + £ or* . 

These equations define a linear velocity field where V 
Is a linear vector function of position# This is termed 
a field of homogeneous deformation. 

2# The Equations of Linear Fields of Motion# 

Keglecting for the time being the translational 
terms of equations (A 1-4), the remaining terms can he 
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written in tensor form. Let the tensor that determines 
the above velocity field be f, and let the conjugate 
of $ be $ c . Then 


U 1-5) 


§- 


fo* 


> 

or, , 

I j 
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o-», ) 


*1 
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CL zx 



aij 

u4j 
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CLlz 



The equations (A 1-4) involving the components of 
V become 


^ + ^-ii X + <in ij + a u 2: 

f A 1-7) . cr= tr 0 + o. 2l * + a lx _ lj + <x n £ 

ur= ar 0+ a. 5l x + a si<j + d ai2: . 

The tensor § can be split into a symmetrical part 
of six terms and an antisyrnnetrical part of three terms. 
The symmetrical part represents a velocity of pure 
extension along the three principal axes, while the 
antisymmetrical part represents a velocity of rotation. 
The symmetrical part is 4(§ + $ e )» and the antisymmetrical 
part is 4($-$c) • 

It can be shown that 

(A 1-8) $= X7 V £ $ c - \/ 7 # ** 

By direct relationship of the above equations, it 
is readily verified that 


** The above developments of sections 1 and 2 are those 
of Frandtl and Tietjens (7). 
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U 1-9) 


<*|| + 4,0 ^(a S i+a,s) j 

^($ + #0=^ %(4n+4ti) a* z. ^(<a> t +a»0 / 

L ^( Au+aii.) a si 


o 4(a lt - an) ^ 

(A 1-10) %($-& ,)- 3*(a,t-a,.) o ^(a te -a»0 

l^/£(a, 4 - a lt ) /i(a ti -a%J o 

Choosing axes so that 


(A 1-11) ^(flw + aa)' /4(<*■!» +a»i)- /i(a 2 5 +a.»t) * o 

then ^($+io) can he written as 


(A 1-12) 


*Z($+ $c) = 



By definition* the axes so chosen are the principal axes. 

Finally, the equations (A 1-4) involving the components 
of V become 


<J-u. + (aS + +CL-45 ~ #2 i)m + 'iCu-ij- 4%i)i: 

{A 1 -13) cr = o; + (a .1 + aA + a& 3 )<j -(a.,' + a«)<j + £(a w - ««)x + £(<*«-a*^ 

W = tui fd» + Oia. + U n ) 2; ~ Cd<i + flzOi + ~^is)x +" &CO-ii- . 


For simplification, new symbols are introduced* 


( 35 ) 




U= (Jo +(a+6+c)x - (e+c)x + D«j + £2 


(A 1-14) 

•j 

(/= Uo + 6A+B+C)ij - (A+C)ij - D X - Fi 
uS-uS a + (A+ £>+ c)£ “ C^ + 8 )2~£X-^y 

where 


A = O x D = - 

(A 1-15) 

B= o', £ = + 

C= uS* F = - 


The equations (A 1-14), involving only nine constants* 
describe the motion of the most general homogeneous 
lineal velocity field. 


\ ■ 
\ 
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j&PPEBDIX II 


Solutions of Linear Differential Equations with 
Constant Coefficients* 

ThiB section applies to the solutions of the 
differential equations of streamlines when the general 
field of motion is considered involving all of the\ 

constants of equation (A 1-14) except translation* The 

\ 

differential equations are \ 


(A 2-1) 


_dx___djj___ 

AxfD«j+£2 B'j-Dx+Fz 


which are equal to 


(A 2-2) 


_ X dx + /<• + -0 di _ ■.\ 

(^A-/iD-n) £) X + (X D+/-1 d~i)F) <-] + (A £ +/J.F+ 


where A #/ u,t> are constant multipliers* 

The expression (A 2-2) will he an exact differential 

\ 

If the numerator is a constant multiple of the differential 


of the denominator} that is 


XA-/x D--$£=pA 
(A 2-3) AD+a 

X £+/x F+ t)c = 

These homogeneous linear equations in the multipliers 
may be solved only if the following determinant 
vanishes t 

A-p -d -£ 

(A 2-4) D B -p -F = o 

£ f C-p 
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This determinant represents a cubic equation that can be 
written as 

(A 2-5) f> 3 -(A+e+c)p*+(A0+AC+6C+D%-£ , '+r ‘){>-(abc + a^+ss*-*- 

S*or each of the roots p, » p 4 , p s of this cubic, the three 
equations for the multipliers may be solved* Let the 
resulting denominators of the differential expression 
(A 2-2) be n , i xi 3 to give 

(A 2-6) d(iyvil 5 ) . 

Then the solutions are 

(A 2-7) - K« . 

The developments of this section to this point are 
by Ford (3), 

A direct solution of the cubic in terms of the 
letter constants is a very laborious algebraic process* 
For an idea of the possible streamline families, an 
indirect solution is much easier by assigning arbitrary 
values to the constants that describe the various types 
of wind fields in the atmosphere* A nomograph by Lipka 
(4-), Alignment Chart for Solution of Quadratic and Cubic 
Equations, can be used for the solutions of the cubic 
where real roots are involved. 

It is found that the solutions of the cubic may 
involve ell real roots* The family of streamlines can 
be found quite readily* The other possibility is that 


co l ) = o . 


•V 
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the solution of the cubic involves one real root and two 
complex conjugate roots. In this case the procedure as 
outlined below must be followed* 

Let f> be the real roots and let 

a + 6-i 

ft* a ' (ri 

where i in this section refers to the square root of 
minus one. Then 

il, = si, (x,y,z) Beal terms only. 

JR.*. ■ [(c+dOx+Ce+fO^ + ^+m)*] 
j\ s = [(c-di.)x + (e-f i)«j +(g-hO^] 

where c,d,e,f»g,h are algebraically complicated functions 
of a and b. 

The solutions of the differential equations will be 
of the form 

(A 2-8) KhIU^“- . 

It can be shown that the third feetor of these equations 
gives the same family of streamlines when used with the 
first factor as doeB the second factor, and hence it 
will not be considered further. 

Taking logarithms of both sides of equation (A 2-8), 
(A 2-9) in-a,= An Ku 4 ^[( r cx+etj + gz)+c (dx+f‘j+hz)] 

But ( p 4 i cjj) = iw A. + t(e + 2 w.tt) . 

Therefore, 
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( A 2-10) 


■jp' inil, = + /£ ^n{j(cx+€y4<3a) z +('dx+f»j + h2.) i ] 


+ i a<*,c Van + h% 

cx+ea+'Jfc 

Separating real and imaginary terms gives 


(A 2-11) 


i>vil, a in K'u + & in ((cx+ey+ga) + (dx+f<j+l'**) 1 ] 


A simultaneous solution of these last equations gives 
the family of streamlines when complex quantities are 
involved. 
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./PPEBDIX III. 


Frontogenesis in Linear FieldB of Motion. 

This section is intended to furnish the intermediate 
steps between equations (9-2) and (10-2). Equation (9-2) 
is 

(A 3-1) f - - ^s' Ws- v) V/f VS j 

where yS is the gradient of any conservative property. 

It can be written as 
(A 3-2) v S = |vsl ( j + yi k) 

where l»m»n are the direction cosines of the vector 
representing the direction of the gradient of the property. 
The equations for linear fields of motion (2-2) can 
be rewritten as follows: 


{A 3-3) 


Y=(Uo + AX+D»j + £2:)c4- (Co + Bvj - Dx + ^i) j 
+ ( llTo + Ci - £x - ^«j) K - 


Substituting this equation into equation (A 3-1)» the 
following steps take places 

VS*0 = |v5|(i '5* + >rv + >\ ^a) 

(vs-v) Y= I*SI[(4 0 X + VV1 + h U*) c 

+ (jL Ux + >w tfy + v\ cq) j 
+ (i Ufx + •tS'y + T) uTt) «J 

-f = -|vS|[(jf + O^+fnUi) 

+ (im Ox + Ojj vn>, Oi ) 

+ ( f n afx + vri-n + -n 1 

Replacing the partial derivatives by their letter constants 
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of equations 

(A 1-15) » the above can be rewritten as 

f= - |vS | \_i-f A + D + in £) 

(A 3-4) 

+ £-ivn D+ vn* B+'n'xn fO 
+ (~j2>v £ - win f 4- n.* - C )J 

or. 

V 

v 3 

\ .V 

(A 3-5) 

■f *-|*S|[fV\ + Vn a B+ n l c] \\ 

or. 

! ' 

(A 3-6) 

•f= |vS|[(e + c)i*+(A*c)w*+(A+B)n 4 - fa+B+c)] 
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